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1. INTRODUCTION 

At present, in majority of the developing countries is known for its trendy, massive, and smart 
commercial complex powered by large number of resources in order to make it ready for production floor. 
One easiest means of automatically controlling, monitoring, and managing the existing resources in 
commercial building is by deploying occupancy sensor [1]. At present, there are various forms of occupancy 
sensor e.g. electromagnetic signal, image sensor, carbon-dioxide sensor, ultrasonic sensor, power meters, etc 
[2], [3]. The methods employed by such sensors are both terminal and non-terminal wise [4]. The base 
function of occupancy sensor is to involuntarily control the switching system of the resources in case of 
unoccupied spaces. At present, there are two dominant techniques for sensing the level of occupancy i.e. by 
using infrared sensor and ultrasonic sensor [5]. Basically, usage of infrared sensors resides in identifying any 
form of change in temperature of the monitored spaces whereas the utilization of ultrasonic sensor is found 
more on acoustics with high frequency in order to perform motion detection [6], [7]. There is also certain 
hybrid level of sensors that uses both of these techniques for enhancing the detection rate with better 
accuracy. However, such sensors are quite expensive in nature. There are certain set of research work which 
has highlighted that positioning of occupancy sensor is extremely important for its successful operation. At 
present, the classification of occupancy sensing process is essentially categorized in three forms (a) based on 
function, (b) based on method, and (c) based on infrastructure. [8]. However, it is normally noted that 
existing occupancy sensing mainly uses sensors where ultrasonic and microwave sensors are much preferred 
owing to its equivalent operation in the form of radar with Doppler shift [9]. Basically, Doppler effect is 
utilized in Doppler radar in order to generate the correct positioning of any object with its specific distance 
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using microwave signal. Usage of Doppler radar is more in aviation industry, meteorology, healthcare, 
ballistic, etc. Maximum usage of Doppler radar is less seen in occupancy sensing model owing to wavelength 
range, longer distances of operation, lack of enhanced frame rate, etc. Even it offers significant amount of 
problems in obtaining less than 1% of error in detection rate. Apart from this, performing system-on-chip 
operation along with compression followed by analysis is still an open challenge with non-compatibility of 
field upgradable. At the same time, there is also a contradiction of theory and practical usage. Occupancy 
sensor was originally meant of resource (energy) consumption but in reality user doesn’t prioritize much on 
energy and emphasize more on their comfort. Majority of the upcoming application calls for instantaneous 
detection of physiological factors of subject where demands of accuracy is more dominant over any other 
issues e.g. resource consumption. We strongly believe that usage of Doppler radar can be truly used for 
enhancing the performance of occupancy sensor. The research area is quite novel as there is no much 
commercially success product on this idea of accuracy for occupancy sensing. Moreover, there is not much 
work being carried out towards exploring the direction of mobility of an object in occupancy sensing using 
Doppler radar. Therefore, we present an analytical model where Doppler radar as well as occupancy sensing 
is integrated with an objective of accomplishing better accuracy. The significant contribution of the proposed 
system is actually to incorporate more granularities in the mobility factor of the subject being traced by 
occupancy sensor using Doppler radar. Section 1.1 discusses about the existing literatures where different 
techniques of both occupancy sensing as well as utilization of Doppler radar are discussed for detection 
schemes used followed by discussion of research problems in Section 1.2 and the proposed research 
methodology as a significant solution in 1.3. Section 2 discusses about analytical model along with algorithm 
description followed by discussion of result analysis in Section 3. Finally, the conclusive remarks are 
provided in Section 4. 


1.1. Background 

This section discusses the existing research contribution towards enhancing the design of occupancy 
sensing. Eedara et al. [10] have presented a technique of sensing occupancy using metadata of standard WiFi 
interface and received signal strength. Study towards the dual-band based input-output based approach of 
sensing subject is carried out by Iyer et al. [11]. Liu et al. [12] have used pyroelectric infrared sensors to 
perform occupancy sensing using standard hidden Markov Model. The work carried out by Li et al. [13] has 
introduced deep sensing technique using the recursive mechanism of a mathematical model. Cooley et al. 
[15] have used a special form of capacitative sensor to perform detection irrespective of any illumination 
condition. George et al. [16] have presented an integrated prototype framework using inductive as well as a 
capacitative approach for designing occupancy sensor. Hossain and Champagne [17] have addressed the 
similar sensing problem with the aid of wideband spectrum considering a case study of cognitive radio. It has 
been seen that there are not much implementation using Doppler radar. On the other hand, the studies 
towards Doppler radar have addressed a different set of problem. Most recently, identification of the target is 
carried out by Cooper et al. [18] using an analytical model. Doppler Radar has also been reported to be used 
in the identification of wall breathing as seen in the work of Aversano et al. [19]. Implementation of Doppler 
radar towards phase-locked loop has been carried out by Mercuri et al. [20] to perform non-invasive 
monitoring of vital signs. Saho et al. [21] have used Doppler Radar for classification of the gaits of a 
different form of subjects in the form of simulation. Wang et al. [22] have presented a scheme for identifying 
the wind component using Doppler radar. Zhao et al. [23] have carried out an implementation of Doppler 
radar to perform detection of physiological characters mainly emphasizing on the low-power transceiver 
design. The problems associated with the correction of the radius in Doppler Radius is addressed in the work 
of Gao et al. [24] where radar of continuous wavelength of high precision is used. Gurbuz et al. [25] have 
used a similar form of waves to perform detection of simple mobile objects using radar signatures. Ritchie 
et al. [26] have performed feature extraction using multi-static waveforms of Doppler radar and was used for 
classifying flying objects. Identification of respiratory patterns using Doppler radar was carried out by Lee 
et al. [27]. The usage of Doppler radar of high frequency has also been reported in the work of Liao et al. 
[28]. The identification of various significant information from the satellites has been carried out by Vega 
et al. [29] where the Doppler radar is also equipped with double frequency as well as polarization. The work 
carried out by Ren et al. [30] has performed identification of heartbeats using demodulation of the complex 
signal for precise detection performance. Majority of the existing approaches using Doppler radar has been 
tested in hardware-based approach with a lesser computational model while the majority of the studies 
focuses on identification of the certain targets. The extracted research issues are briefly outlined in next 
section. 
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1.2. Research Problem 
The significant research problems are as follows: 
a. Existing techniques are much less on using Doppler radar for assisting in occupancy sensing 
modeling and implementation. 
b. Less computational approach and more hardware-based approach towards Doppler radar usage 
render reduced scope in the presence of uncertainty and dynamic modeling. 
c. Existing framework is capable of identifying the presence of certain objects; however, its 
granularity towards investigating the physiological characteristics is completely missing. 
d. A robust mathematical model for integrating occupancy sensing and usage of Doppler radar is 
few to find in existing approaches. 
Therefore, the problem statement of the proposed study can be stated as “Designing an approach for 
offering better granularity in the mobility behavior of a moving object for enhanced occupancy sensing by 
harnessing the potential of Doppler Radar." 


1.3. Proposed Solution 

The prime purpose of the proposed system is to perform an effective modeling of applying Doppler 
Radar for better occupational sensing. The contributions in the present work are a) inclusion of dynamic and 
uncertain mobility of subject b) incorporation of granularity of mobility by incorporating directionality 
factor. By the term granularity, it will mean that the proposed system not only perform detection of the 
presence of any subject within the specified range of Doppler Radar but also computes its rate of specific 
movement that is quite periodic. For simpler model construction, we consider to compute the rate of 
respiration for the subject as well as its direction towards and away from the reference point (where the radar 


resides). The flow of the proposed method is shown in Figure 1. 
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Figure 1. Proposed flow of Methodology 


According to the proposed system (Figure 1), we first simulate the Doppler radar in order to perform 
more granular investigation of the biological signal by reading its input signal followed by generating the 
input signal performed by the filtering operation as well as implicating the framing & windowing operation 
followed by the accomplishment of the local maxima and FFT operation. The prime logic of proposed study 
is that there is a significant amount of shifting of frequency over the spectrum of input signal whenever a 
subject performs any motion. Depending on the velocity of the mobility of the subject, the direction, as well 
as score of frequency, can be computed. The proposed system, therefore, emphasize on the received signal 
spectrum. We also apply a simple modulation-based framework that considers uncertain mobility issues to 
trace the spectrum with shifted frequency. Hence, in a nutshell, the study will finally aim towards observing 
the shifted frequency of the spectrum to compute the rate of any action related to the biological signal. An 
analytical model is designed for this purpose where simpler mathematical modeling is applied for 
constructing the identification of the subject as well as computation of the rate of respiration under the very 
uncertain condition of mobility. The next section discusses the analytical model. 
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2. ANALYTICAL MODEL 

This section discusses the analytical model being implemented towards designing a novel 
occupation sensing system for uncertain and undynamic subject mobility using Doppler radar. The presented 
analytical model makes use of the Doppler shift to track the presence of any forms of biological signals of 
any object. The applicability of Doppler radar is quite high in this as according to theory, we can say that 
there is a periodical mobility associated with an object with the same frequency. However, the object might 
have a phase shift with time-varying nature, and therefore, it is much suitable for detection of any subject 
with the biological signal in a non-invasive manner. A simple example to understand this analytical 
implementation is the process of respiration that is characterized by periodic displacement of the chest. 
Moreover, such mobility factors will introduce a highly varying phase shift better proportion to the chest 
displacement considering signal wavelength. The proposed system applies modulation scheme of Binary 
Phase Shift Keying (BPSK) although other forms of modulation could also experiment. The algorithm 
develops a channel with MIMO form to construct channel object. The prime emphasize was offered in 
developing a cost-effective implementation of identification process of biological signals of any inhabitants 
of the specific area. Another uniqueness of proposed algorithm is that it doesn't offer any form of recursive 
operation. Thus it leads to faster response time as well as lower computational complexity in the form of 
storage complexity. All the computation is carried out at run-time, and there are no residual signals being left 
behind in memory. The complete modeling is carried out using six discrete operations, i.e., performing a 
simulation of Doppler radar, generating an input signal, performing filtering operation, implicating framing 
and windowing operation, obtaining local maxima, and finally performing Fast Fourier Transform (FFT) 
operation. Following are the brief outline of the proposed analytical model. 

a. Performing Simulation of Doppler Radars 

This algorithm is essentially meant for simulating the Doppler radar that is deployed for detecting 

the presence of any new subject by their biological signal. The algorithm takes the input of S 

(stream), M (order of modulation), N, (transmitter antenna), N, (receiver antenna), fy (maximum 

Doppler shift for all path), t, (input sample period), t (path delay), Ppa (mean path gain), o (antenna 

coefficient), Nsam ¢ (Samples per frame), and N,amp (Samples per frame) that after processing gives the 

output of h, (Doppler spectrum). The steps included in the algorithm are: 


Algorithm for Simulating Doppler Radars 
Input: S, M, Na No, fa, ts, T, Poa, O, Nsam f Nsamp 
Output: h, 

Start 

1. init S, M 

2. h=mc(N,, N;, ts, fa, T, Poa) 

3. h.TxCM>a(0) 

4. For iF=1:N frame 

Sigin=Mod(Amodem; arb(S, M), Nsam f N 
idx=( 1 :Nsam_ +F- 1 )*Nsamp 

filter(h, Sigin) 

For i=1: N, 

9. [yi y2, Y3, (idx, it)]J=h.PgainGit) 

10. End 

11. End 

12. Show PSD(hs, ‘Nsamp/5”) 

End 


POSTON 


The first step of the algorithm is to develop a set of local random number stream using Modified 
subtract with borrow generator with a seed for generating the local stream S (Line-1). The algorithm 
considers 2™ order of modulation M to create an object of BPSK modulation hyogem(Line-5). The 
input sample period ts is computed as follows: 


1 
pce C= (1) 


s N,. 
In the above expression, the input bit rate is obtained by the product of input symbol rate R sym with 
the logarithmic function of modulation order M. The variable N,, is a factor of oversampling. The 
next part of the algorithm implementation is associated with the channel modeling by considering a 
number of the factors as shown in Line-2. The essential factors are path delay t, mean path gains (in 
dB) Pea, maximum Doppler shift allocated for all paths fa, some transmitting / receiving antenna 
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N/N,, correlation coefficient that is equivalent to the antenna correlation o. The channel object h is 
created (Line-2), and we also derive a transmit correlation matrix shown in Line-3 using non- 
symmetric Toeplitz matrix over the argument of antenna correlation. The study computes all the 
number of the frames Nérame Obtained by dividing a total number of channel samples N samp by some 
samples per frame Ngamp ¢- A loop is created for all the number of frames (Line-4), where the prime 
task is to compute input signal Sig;,, index idx, and output. The input signal is formed by applying 
modulation function on modulation object hmodem, an arbitrary value of signal and modulation 
order, sample frequency, and a number of the transmitter (Line-5). The index is formed by adding 
all the individual sample frequency with the product of remaining frames with sample frequencies 
(Line-6). The output is formed by applying filtering operation on channel object h and input signal 
Sigin (Line-7). Another nested loop (Line-8) is created for all the number of transmitters where the 
relative information of path gains Pain concerning channel object h is computed respectively for yj, 
y2, and y3 matrix (Line-9). Finally, a power spectral density psd is obtained (Line-12) to obtain the 
simulation of Doppler radar. 


b. Generating Input Signal 
The input signal is formed by considering synthetic biological signal data and representing it in the 
form of time (x-axis) and amplitude (y-axis). 


c. Performing Filtering Operation 
To perform filtering operation, the algorithm considers inclusion of a new variable 6 that stands for 
preemphasis coefficient. The filtered outcome F is obtained as follows: 


F =((1-6],x,;.) (2) 


In the above equation, the function ¢ represents a filtering operation on Xpio (vector derived from the 
input signal) with the filter represented by [1- 6] and 1 in order to obtained filtered data F. The 
preemphasis coefficient 6 is used for maximizing the measure of certain frequency of higher order 
for enhancing the cumulative signal quality. It is also responsible for controlling any form of 
possible effects of signal saturation, distortion, or attenuation, etc. 


d. Implicating Framing and Windowing Operation 

As the proposed system considers biological signals to be non-stationary therefore for stabilized 
computation, framing and windowing operation is highly essential. This algorithm takes the input of 
ô (preemphasis coefficient), T, (frame duration for analysis), T, (frame shift for analysis), f, 
(Sampling frequency) which after processing yields an output of E (envelop). The steps involved in 
the algorithm are as follows: 


Algorithm for Framing & Windowing 
Input: ô, T,,, Ts, fs 

Output: E 

Start 

1. init ô, T,, To fo 

2. F = (Ul = ô], Xpio) 

3. [Nw Ns] > 1E-w, where y=3*Ts*fs 

4. Frames > 8, , panl F, N, N,) 
5. E=[tmax tmin] DBF, c) 

End 


w? 


The initial steps of this algorithm are same as that of performing filter operation in the prior stage 
(Line-1 to Line-2). The algorithm computes both durations as well as shifts of the given frame 
(Line-3) and then applies a function gẹ frame tO perform framing and windowing operation (Line-4). 
This step is further followed by the computation of peak detection for facilitating windowing 
operation concerning maximum fmax and minimum fmin value of the envelope E (Line-5) considering 
the permissible envelop coefficient of c. 


e. Obtaining Local Maxima 

This process is responsible for computing the local maxima of the given biological signals of the 
Dopler Radar. The primary intention of this step is to obtain the interpolated version of an envelope 
obtained in the prior step. The complete operation of this implementation is similar to Line-1 to 
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Line-4 of an algorithm for framing and windowing with an additional step of applying a function of 
local maxima on obtained frames of Line-4 to correctly obtain upper and lower envelop. 


f. Fast Fourier Transform (FFT) Operation 

This is the last stage of the proposed implementation where FFT operation is carried out on frame F. 
Accomplishment of this step confirmed the form of biological signals. 

We use this model to generate the transmitted signal, that is normally unmodulated can be 
represented empirically as: 


T signa Cos (2nft+ (7/180) *30)) (3) 


The next phase of the implementation is to construct a received signal. We formulate an algorithm 
that takes the input of c (velocity of light), fı (frequency), a (heart mobility factor), b (respiratory 
mobility factor), ¢,(time), v (velocity of object), dini (initial distance) in order to yield an output of im 
(integrated mobility), R (received signal). 


Algorithm for obtaining the received biosignal concerning Mobility 
Input: c, fi, a, b, ti, v, dini 

Output: im, R 

Start 

1. init dini, V, tı, a, b, fi, € 

2.im= dini+ vti tat b 

3. R=cos((2zft)-((4ad_1)/(c/f,)+ph) 

End 


The complete formulation of the biosignals is carried out by ignoring the variation in amplitude. At 
the same time, we can signify the term vt, as a positive or negative sign for showing the forward and 
reverse the direction of mobility of the subject. The next part of the study is to formulate a spectrum 
for this which is generated by considering a periodic mobility due to heart beat (a) as well as 
respiration (b) in the form of [a b]=sin(2zf,t). According to our formulation, the routine mobility of 
the body is due to the heartbeat and respiration. This method is then continued by computing motion 
modulation where we consider samples per second to compute the wave in the form of: 


x=o(fft(0.5.cos(2zF,.)\+cos(27F;;.t))) (4) 


The function œ represent FFT shift operation that leads to the generation of the spectrum concerning 
frequency and the absolute value of x/N, where N represents the number of the sample concerning 
time. This process offers an empirical response to mobility in forwarding direction concerning 
biological signal. The similar process is applied for computing mobility in a backward direction with 
(t-to) being replaced in the position of t in the above mentioned equation. We also compute the 
actuator result for both before and after shift concerning DFT value and normalized frequency value. 
According to our mobility model, the frequency spectrum will always be inclined towards +ve axis 
if the object is found to be moving towards Doppler radar by 2v/wavelength, while it takes the 
similar empirical value for reverse movement direction as it moves to the -ve axis. This formulation 
of directionality is carried out to resist any forms of uncertainty in the mobility model. The 
advantage of this model is that if the maximum peak is studied, it can give the true picture of its 
shifting behavior towards positive/negative axis which is directly a representation of a minor 
movement even like respiration that can be easily identified using the present model. At the same 
time, the system minimizes the dependencies on multiple radar system as complete information of 
mobility of an object can be tracked by the proposed system concerning its respective direction. 


3. RESULT ANALYSIS 

The numerical values considered to carry out this analysis are: preemphasis coefficient 6 of 0.01, 
modulation order M=2, input of symbol rate 10000, oversampling factor N,.=4, maximum Doppler shift 
f,=0.5, 2 transmit antenna, and 1 receive antenna, average path gain of 0 dB, -15 dB, and -20 dB, correlation 
coefficient o=0.7, number of samples per frame Nsamp =1000, total channel samples=1.5x10°. The 
corresponding outcomes are shown in Figure 2. 
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Figure 2. Graphical Outcome of Proposed System 


25 


The graphical outcome shown in Figure 2 shows an effective identification of the biological signal 
that can effectively trace any form of dynamic or random mobility subjects. The complete computation time 
is found to be 0.024487 seconds in core-i3 windows machine. The complete analysis is carried out using 
synthetic respiratory data. The streamlined pattern of the curve in Figure 2(d) to Figure 2(f) shows that 
proposed system can provide significant assists in extracting the rate of respiration with a higher probability 
of lower resource dependencies. The next part of the experiment for mobility considers frequency as 10° Hz, 
the speed of light as 3*10°mps, samples per second as 100, F.=5 Hz, and F,,=0.1 Hz, time to=1 sec. Figure 3 
highlights that proposed system is completely capable of analyzing forward motion (Figure 3(a)) and 
backward motion (Figure 3(b)). Similarly, a better visibility of the spectrum as an impact of shift operation 
can be seen in Figure 3(c) and Figure 3(d). Owing to the non-recursive formulation of the complete technique 
of occupancy sensing in the proposed system, there is a significant scope of better compatibility in low- 
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powered and low-cost hardware design in future. We have also performed computation concerning other 
modulation scheme as well as by changing other simulation variables. However, in the majority of the cases 
there is a uniformity in the detection performance, and thereby it exhibits the better reliability in its outcome 
with an accuracy of 97.85%. Hence, the proposed system offers a cost-effective approach to model 
occupation sensing concerning direction. 
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Figure 3 Analysis of Mobility Patterns 


4. CONCLUSION 

Usage of occupancy sensing, as well as Doppler radar, is not novel from an individual viewpoint, but 
the integrated usage of both is very less reported in existing literature. Although there is an existing system to 
perform detection of a moving object, they miss granularity in modeling the mobility behavior. Therefore, the 
proposed system introduces an analytical model that is used for computing the periodic physiological factor of 
a mobile subject under the monitoring range of the radar. An explicit modulation technique is applied to study 
directionality of object movement using the positive and negative axis of frequency. The study outcome shows 
that proposed system offers better resolution of the waveforms for enriched investigation of the dynamic 
mobility pattern of any object. 
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